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ABSTRACT 

We present results from a study of the orbits of eclipsing binary stars (EBs) 
in the Magellanic Clouds. The samples comprise 4510 EBs found in the Large 
Magellanic Cloud (LMC) by the MACHO project, 2474 LMC EBs found by the 
OGLE-II project (of which 1182 are also in the MACHO sample), 1380 in the 
Small Magellanic Cloud (SMC) found by the MACHO project, and 1317 SMC 
EBs found by the OGLE-II project (of which 677 are also in the MACHO sample); 
we also consider the EROS sample of 79 EBs in the bar of the LMC. Statistics 
of the phase differences between primary and secondary minima allow us to infer 
the statistics of orbital eccentricities within these samples. We confirm the well- 
known absence of eccentric orbit in close binary stars. We also find evidence for 
rapid circularization in longer period systems when one member evolves beyond 
the main sequence, as also found by previous studies. 
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Introduction 



The two components of a binary star system raise tides in each other. These tides will, 
in general, cause dissipation in the stars, leading to the phenomenon of tidal lag. These lags 
exert torques which lead to exchange of angular momentum between the orbit and the spins of 
the stars. The lowest energy configuration, for a given total angular momentum, is a circular 
orbit with both spins aligned and synchronous with the orbit. All other configurations evolve 
towards the synchronous, circular state. Whether or not this is achieved is determined by 
the rate of dissipation in the tides. 

Th e study of tides raised in celestial bodies by mutual attraction dates back at least to 



Darwin! (I1879I ). Darwin considered the tides raised on the Earth, modeled as a homogeneous 
and deformable viscous body, by a point mass Moon in the weak friction limit; hence his 
model is also known as weak friction model. The model assumes that viscous dissipation 
causes a delay in the onset of the tide by a constant amount r so at time t the shape of the 
Earth is the one that it would have been at time t — r in the absence of dissipation. The 
axis of the tidal bulge is not therefore aligned with the line of the centers of the two bodies 
but lags it by a constant amount, resulting in a torque that tends to align the two bodies. 
Modern discussions of the weak friction model, with empha sis o n clo s e bin ary stars rather 
than planet-satellite systems, are given by lAlexanderl (119731 ) and iHutl (11981I ). 

proposed that the coupling of the tidal flow with turbulent flows 



in 



the envelope of a late type star is chiefly responsible for orbital circularization and syn- 
chronization in th ese stars; evidence of cir cularization due to tidal interactions in late type 
giants is fou nd by iLucy fc Sweeney! (Il97ll ). The t heory has been m o re rec ently revised by 
Zahnl (119891 ) and compared against observations by IZahn fc Bouchetl (119891 ) who found that 
for late type stars most of the circularization occurs during the Hayashi p hase of pre main 
sequence evolution, as previously suggested by lMayor fc Mermilliodl (11984 ). 



The circularization models cited above consider only the equilibrium tide which arises if 
the star is at all times in hydrostatic equilibrium. If, however, the orbit is not circular or the 
rotation is not synchronous the star is subjected to a time-varying gravitational potential 
which excites oscillations, giving rise to a dynamical tide superi mposed on the e quilibrium 
tide. Forced os ci llation s in b inary stars were first considered by ICowlingl (Il94ll ) and then 



by I Zahnl (Il970l ); I Zahnl (119751 ) considered the damping of the dynamical tide by radiative 
dissip ation in the radiative envelope as a possib le circularization mechanism for early type 
stars; iGiuricin. Mardirossian. fc Mezzettil (119841 ) studied a sample of ~ 200 early type binary 
stars and showed th e theory to be compatible with the data presented. These mechanisms 
are all described in I Zahnl (119771 ). where time scales of circularization and synchronization 
are derived. 
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A circularization mechanism in which energy dissipation is due to large scale hydro- 
dynamical flows resul ting from the deformation of the star (m eridional circulation) has been 



pr oposed by Tassoul (Tassoulj 



9871; iTassoul &: Tassoullll995l ). This model has been applied 



by lClaret. Gimenez fc Cunhal (119951 ) to a homogene ous sampl e of 45 eclipsing binary stars 



(EBs) with accurate parameter determination from lAndersenl (119911 ): they find a satisfac- 
tory agreement with the observations and also find that circularization is still taking place 
during the main sequence for early-type EBs, a finding consistent with the results of this 



paper. I Claret fc Cunhal (119971 ) then apply b oth the turbulen t dissipation mechanism and the 
radiative damping mechanism to the same (lAndersenl 1 1 9 9 ll ) data set. They find that within 
uncertainties, these formalisms seem able to explain the observed eccentricity distribution, 
although with a few exceptions. 

Since tidal forces decrease and the period increases with separation, in a sample of 
coeval binary stars such as those in star clusters, binaries with a longer period should have 
orbits with a range of eccentricities since tidal forces have not circularized them yet, whereas 
binaries with a shorter period should all have circular orbits since tidal forces would have 
been more effective in circularizing them. Therefore clusters should show a transition between 
binaries with shorter period and circular orbits and binaries with longer periods and eccentric 
orbits; the values of this transition period varies with the age of the cluster. Binary stars in 
cluster are widely studied: for example as a part of the WIYN Open Cluster Study (WOCS) 
0, a systematic search of EBs i n Milky Way open clusters for testing orbital circularizatio n 
theories is currently u nder way (IMathieu. Meibom fc Dolanll2004l ; iMeibom &: Mathieull2005l ) . 
Latham et all (l2002h present spec troscopic orbital solutions for 171 single-lined stars from 



their catalog (jCarney et al.l Il994l ) of 1464 stars selected for high proper motion and find 
that for the metal poor, high velocity halo binary stars in their sample the transition from 
circular to eccentric orbits occurs at ~ 20 days. 

Binary systems in which the component stars eclipse each other are a powerful tool for 
the study of circularization theories since in this case both stellar parameters like radius, 
mass, and temperature, and orbital parameters like eccentricity, angle of inclination, and 
longitude of periastron, can be determined to the accuracy necessary to test different theories. 
To determine these parameters fully one needs to supplement the photometric data with high 
quality spectra, which can be used to determine orbital velocities. 

Eclipsing binary stars have been found in large numbers by astronomical surveys study- 



1 http : //www . astro .uf 1 . edu/~{}ata/wocs/ 
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ing gravitational microlensing: MACHO0, OGLFjfl and EROS0. In particular a sample of 
4634 E Bs in the Large Ma gellanic Cloud (LMC) and of 1509 EBs in the Small Magellanic 
Cloud (jFaccioli et al.l 120071 ) have just been p ublished by the MACHO collaboration. Previ- 
ously a ca talogue of 2580 EBs in th e LMC JWvrzvkowski et al.l [20031 ) and of 1351 EBs in 
the SMC ( iWyrzykowski et al.l 120041 ) were published by the OGLE collabora tion, the EBs 



were selected from their catalogue of variable stars in the Magellanic Clouds (jZebruh et al. 



2001 |) compiled from observat i ons t aken during the second part of the proje ct (OGLE-II: 
Udalski. Kubiak fc Szymahskil 119971 ). A cat alogue of 611 EBs in the LMC ( lAlcock et al. 



1997 ). and of 79 EBs in the bar of the LMC ( Grison et al.|[l~995"l ) were also published by the 
MACHO and EROS collaboration respectively. This earlier MACHO sample is contained in 
the new one. 

We aim to take advantage of the sheer size of our samples to infer statistically valid 
conclusions on orbital circularization, even though our data do not have high photometric 
precision and there are no spectroscopic data. We exploit a simple idea: the phase difference 
between the primary and secondary minima in an EB light curve, 0i — 02, may be effectively 
estimated for the light curves in our samples. This object is related to t he orbital elements 
by th e following equation, which sets a lower bound to the eccentricity (jKallrath fc Milone 
1999h : 



- 02 1 - 



-e cosu I 1 + 



7T 



sin 2 % 



(1) 



where 0i — 02 is the phase difference between primary and secondary minima, e is the 
eccentricity, u is the argument of periastron and i is the orbital inclination. For eclipsing 
systems, especially the wider systems where non zero eccentricities are found, we know that 
sini ~ 1. The important degeneracy is between uj and e. It is possi ble in principle to 



mode l high signal-to-noise photometric da ta to break this degeneracy (IWilson fc Wyithe 



200ll ; IWilson fc Wvithel 120021 ; iDevorl 120051 ) but we do not attempt this here. We will show 
that 0i — 02 alone may be used for our purpose. Although Eq. [1] shows that for an eccentric 
orbit 0i — 02 can be either < 0.5 or > 0.5, depending on cosu, we will always adopt the 
convention 0i — 02 > 0.5 since only the deviation of 0i — 02 from 0.5 is relevant for the 
detection of eccentricity. 

The paper is organized as follows: Section [2] describes both the MACHO and the OGLE- 
II samples for both Clouds; Section [3] describes the fits to the MACHO data and its validity 



^htt pT//www. macho .monaster . ca/| 

3 http : //sirius . astrouw . edu.pl/~{}ogle/ 
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(for the OGLE-II data the relevant information is provided by the authors); Section H] reports 
our results for both Clouds both from the MACHO and the OGLE-II samples; Section [5] 
discusses the significance of our results and EROS results are briefly considered in Subsection 
15.21 Finally Section states our conclusions. 

The data presented in this paper can be accessed on line at the Astrophysical Journal 
website^] and are mirrored at the Harvard University Initiative in Innovative Computing 
(IIC) /Time Series CenteiB 

From now on we will always use the term unfolded light curve to indicate a set of time 
ordered observations and will reserve the term light curve to indicate a set of time ordered 
observations folded around a period, omitting for brevity the adjectives "folded" or "phased" . 



2. Data sets 

2.1. MACHO data 

The MACHO Project was an astronomical survey whose primary aim was to detect 
gravitational microlensing events of background stars by compact objects in the halo of 
the Milky Way. The background stars were located in the LMC, SMC and th e bulge of the 



Milky Way; more details on the detection of microlensing events can be found in lAlcock et al. 



(120001 ) and references therein. 



Observations were carried out from July 1992 to December 1999 with the dedicated 
1.27m telescope of Mount Stromlo, Australia, using a 2 x 2 mosaic of 2048 x 2048 CCD 
in two band passes simultaneously. These are called MACHO "blue", hereafter indicated 
with Vmacho, with a bandpass of 440 — 590nm and MACHO "red", hereafter indicated with 
-Rmacho; with a bandpass of 590 — 780n m; transforma t ions t o standard Johnson V and 



Cousins R bands are described in detail in lAlcock et al.l (119991 ); the magnitudes quoted in 



this paper have been obtained by using the following transformation for the LMC: 

V — ^macho + 24.22 mag — 0.18(Vmacho — -Rmacho) 

R = -Rmacho + 23.98 mag + 0.18(Vmacho — -Rmacho)- (2) 
and the following one for the SMC: 

V = Vmacho + 24.97 mag - 0.18(Vmacho - -Rmacho) 



'http : / / www . j ournals . uchicago . edu/Ap J/ 

e http : //timemachine . iic .harvard. edu/f accioli/CircularizationTables/ 
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R 



-Rmacho + 24.73 mag + 0.18(Vmacho — -Rmacho)- 



(3) 



The zero point in Eq. El is different from the one in Eq. [2] because of the different exposure 
times for LMC and SMC (Ulcock et al.lll99flh . 



There are several hundreds of observations in both band passes for most EBs; the 
central fields of the LMC were observed more frequently than the periphery and there are 
on average fewer observations in the red band because one half of one of the red CCDs was 
out of commission during part of the project. 



For each source in the database found t o be variable (jCook et al.lll995l ) a p eriod was 
found using the Supersmoother algorithm (iReimannl 1 19941 . first published by iFriedman 
(1l984j )). The algorithm folds the unfolded light curve around different trial periods and 
selects the one that gives the smoothest light curve. Periods were found separately for the 
red and blue unfolded light curve and they usually agree to hig h accuracy (0.03% o n average); 
unfolded light curves were then folded around these periods (jFaccioli et al.l 120071 ). 



2.2. OGLE-II data 



The OGLE-II d ata w e considered comprise the LMC sample of 2580 EBs described in 
Wyrzykowski et al.l (120031 ) of wh ich 1182 were also presen t in our sample, and the SMC 
sample of 1351 EBs described in IWyrzykowski et al.l (12004 ) of which 677 were also present 
in our sample. The available data include the phase of secondary minimum <ft sec (the phase 
of primary minimum is set to 0), the period, and B, V, I, and Jdia magnitudes at maximum 
l ight. Here B, V, and / refer to magnitudes obtained using standard PSF fitting photometry 
(jSzymariskil 120051 ) and Jdta are I band magnitudes computed via Differen ce Image Analysis 
photometry (DIA: IZebruh. Soszyhski. fc Wozniakl|200ll ; ISzymahskil 120051 ) . We used both / 
and /dia magnitudes in our analysis. The phase of the secondary sec is given up to two 
decimal places. 



3. Light curve fitting 

To find the phases of minima for the MACHO samples the light curves have been fitted 
to a sum of sines and cosines: 

m(4>) = ^j{cLj cosujj(f) + bj sin ujj(p) (4) 

where m(<j>) are the instrumental MACHO magnitudes, Vmacho and -Rmacho in which the 
fits have been carried out and <fi is the orbital phase for the light curves. We preferred 
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instrumental MACHO magnitudes because in many cases observations in one of the MA- 
CHO bands are invalid, thus forcing us to exclude these points in the fit for both bands if 
standard magnitudes ha d been used . The "frequencies" Uj have been found with the Lomb 
Periodogram technique ( Lomb 1976 ; Scargle 1982 ; Press et"aD 1992 ); this method is not a 
Fourier decomposition because the frequencies are not in a harmonic series; furthermore, the 
data are not uniformly spaced in phase. Before fitting, outlying points in the light curves 
have been eliminated using moving windows with roughly ~ A^tot/50 (where ./Vtot is the 
number of points in the light curve) points, calculating their mean and standard deviation 
and excluding from the fit points that are > 2 standard deviations away from the mean in 
each window. Figure [1] shows some examples of sinusoidal fits and data on these objects are 
shown in Table [TJ The phasing was performed using the red and blue periods for the red and 
blue bands respectively. The number of frequencies calculated by the Lomb Periodogram 
varies from object to object, but in most cases is in the range 10 — 100. 

One problem of sinusoidal fit is that it tends to underestimate the depth of the minima 
for deep, widely separated eclipses; in general this fit works better for tidally distorted sys- 
tems and less well for widely detached systems. Therefore, after performing the sinusoidal 
fit, the light curves were also fitted to a sum of two Gaussians, which often better represent 
detached systems. We used the positions of the minima determined via sinusoidal fit as a 
startin g point for the G aussians which was performed via the Levenberg-Marquardt algo- 
rithm (IPress et al.lll992l ). The best of the two fits from either sinusoidal or Gaussian fit for 
each light curve in each band was then selected. Figure [2] shows some examples of Gaussian 
fits; data on these objects are shown in Table [TJ 

It is necessary to estimate the errors in our determinations of the quantity (pi — 02- 
We did this as follows for each of our systems. First, w e created a synthetic light curve 
using the EBOP code (jEtzell Il98ll ; iPopper fc Etzellll98ll ). which implements the model by 
Nelson fc Davi] Jl972h with some modifications and we followed the prescriptions described 



in 



Alcock et al. J 997h . his synthetic light curve was used only to test our procedure for 
estimating (pi — (p2- and was sampled at the same phases (pi as the original light curve; noise 
was added to mimic a MACHO light curve. The synthetic noise was determined from the 
EBOP fit to the MACHO light curve as follows. First the residuals Oj — Cj were determined, 
where Oi is the observed magnitude at phase (pi an d Cj the value of the fitted magnitude 
at that phase. These residuals were divided by the MACHO estimated photometric errors 
<jj to create scaled residuals Si = 0i ~ Cl . The cumulative distribution of the Sj was sampled 
randomly and used to create the synthetic noise at phase (pi, Ni — SjO~i. The new, synthetic, 
noisy light curve was then Of — Ci + N{. We determined the quantity (pi — (p2 for this light 
curve, and the entire procedure was repeated ~ 30 times in each band. We excluded the 
cases in which (pi — (p 2 > 0.9 (as we did in the real fits) which we took as an indication that 
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Fig. 1. — Examples of light curve fitting with trigonometric functions used for minimum de- 
terminations for some LMC EBs in the MACHO sample. Triangular points indicate outlying 
points excluded from the fits. The light curves are shown in order of increasing period. 
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Fig. 2. — Example of light curve fitting with Gaussians used for minimum determinations 
for some EBs in the LMC. Triangular points indicate outlying points eliminated from the 
fit. The light curves are shown in order of increasing period. 
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the fitting procedure had failed. The RMS of the deviations from 0i — 02 determined from 
the MACHO data is a reasonable estimate of the error in this procedure. Figure [3] shows 
the histograms of the errors thus defined for both bands; most errors are smaller than 0.01 
making our method for detecting eccentricities very robust. We then studied the effects of 
period uncertainties on our method. A small period uncertainty should manifest itself as 
a constant shift in the phases of a light curve: therefore we selected a random subsample 
of ~ 50 EBs, created simulated light curves in the same way but in addition we added a 
constant random phase shift to them and repeated the procedure. The results are shown 
in Figure H] which shows that our method is also robust with respect to uncertainties in the 
period determination. 

In the few cases where our procedure gives a large error estimate the reason is mostly 
the presence of one or two shallow minima and/or noise in the light curve. In the presence 
of a very shallow eclipse a small perturbation of the data can lead to a large difference in 
the determination of one or both phases of the minima and therefore a large variance used 
to estimate the error in phase difference; this is true also for noisy data. 



4. Results 
4.1. LMC results from MACHO 

For the LMC our fitting procedure gave acceptable results in at least one band for 



4510 EBs out of the 4634 EBs that make up the LMC sample described in iFaccioli et al. 



( 120071 ). Figures [5] and [6] report the Color Magnitude Diagram and the Color Period Diagram 
respectively for the LMC EBs in our sample in different panels according to their phase 
difference. The explanation of the panels is given in Table [2], which also reports the minimum 
eccentricity corresponding to a given phase difference, obtained by assuming cos a; = 1 in 
Eq. [TJ For the purpose of this discussion we describe EBs with V — R < 0.2 mag as Main 
Sequence systems and EBs with V—R > 0.2 mag as evolved systems; we also describe EBs 
with P > 20d as "long period" EBs; the breakup of MACHO LMC and SMC EBs and 
OGLE-II LMC EBs used in our analysis is given by Table [31 

It is evident from the Color Magnitude Diagram (Figure [5]) that evolved systems have 
102 — 0i | < 0.6 in the vast majority of cases; indeed there is a cutoff in phase difference at 
V — R ~ 0.2 mag. The Color Period Diagram shown in Figure [6] shows the cutoff in V — R 
too, and also shows that elliptical orbits are concentrated in the 1.5d < P < 20d range. 
This is not surprising in view of Kepler's third law P oc a 3 / 2 , where a is the semi- major axis: 
in a system with long period and hence widely separated stars, tidal forces are less effective 



Table 1. Basic data for LMC EBs with light curves shown in Figures [T] and [2j The light 
curves are shown in order of increasing period. 



MACHO ID 


RA(J2000) 


DEC(J2000) 


P(d) 


ya. 


V-R a 


01 - 02 


Type of fit 


17.2586.264 


04:56:08.462 


-69:56:34.48 


1.56 


16.50 


-0.06 


0.50 


Sinusoidal 


10.4035.145 


05:05:02.233 


-70:06:13.27 


2.53 


17.24 


-0.05 


0.55 


Sinusoidal 


11.8867.249 


05:35:14.688 


-70:40:19.12 


3.14 


17.72 


0.01 


0.50 


Gaussian 


1.3804.164 


05:03:36.536 


-69:23:32.27 


4.19 


17.20 


-0.02 


0.50 


Sinusoidal 


79.4899.724 


05:10:21.705 


-69:01:01.48 


5.00 


18.05 


0.03 


0.62 


Gaussian 


82.8285.11 


05:31:04.061 


-69:09:20.81 


11.61 


16.85 


-0.04 


0.50 


Gaussian 


78.5618.83 


05:15:19.305 


-69:26:39.97 


15.96 


14.35 


-0.04 


0.61 


Sinusoidal 


22.4989.1285 


05:11:31.896 


-71:01:45.95 


107.24 


18.28 


0.55 


0.74 


Gaussian 



a Values are quoted to the hundredths of magnitude, typical of MACHO observational uncer- 
tainties. 



Table 2. Explanation of panel labels for Figures El EH H El EH HUHTU and [TJ 



Panel label 


Phase difference 


Minimum eccentricity 


(a) 


102 - 0l| < 0.51 


e > 


(b) 


0.51 < |02 - 0i | < 0.6 


e > 0.016 


(c) 


0.6 < |02- 0i | <0.65 


e > 0.16 


(d) 


|02 - 0i | > 0.65 


e > 0.24 



Table 3. Summary of the sample features. 



Galaxy 


Sample 


Total 


MS a 


Evolved b 


Long Period 


Long Period MS a 


Long Period Evolved b 


LMC 


MACHO 


4510 


3667(81%) 


843(19%) 


349 


21(6%) 


328(94%) 


LMC 


OGLE-II 


2474 


1744(70%) 


730(30%) 


216 


31(14.3%) 


185(85.7%) 


SMC 


MACHO 


1380 


1293(94%) 


87(6%) 


66 


19(29%) 


47(71%) 


SMC 


OGLE-II 


1317 


1091(83%) 


226(17%) 


165 


27(16.4%) 


138(83.6%) 



a Main Sequence: defined as V — R < 0.2 mag for the MACHO samples and as V — Idia < 0.4 mag for the 
OGLE-II samples. 

b Defmed as V — R > 0.2 mag for the MACHO samples and as V — Idia > 0.4 mag for the OGLE-II samples. 
°Defmed as P > 20d both for the MACHO and the OGLE-II samples. 
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Fig. 3. — Histograms of the errors in the phase differences of the minima from the Monte 
Carlo simulation. 
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Fig. 4. — Histograms of the errors in the phase differences of the minima from the Monte 
Carlo simulation with a random shift added to the phases. 
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Fig. 5. — Color-Magnitude Diagram for 4510 EBs in the LMC MACHO sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i] < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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Fig. 6. — Color-Period Diagram for 4510 EBs in the LMC MACHO sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i] < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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at orbital circularization and orbits remain elliptical longer. Conversely in systems with 
shorter period, tidal forces can achieve orbital circularization already on the Main Sequence. 
The diagram shows a strong Color-Period relation for evolved, long period systems. Th ese 
systems are mostly ellipsoidal, with continuously varying light curves (jFaccioli et al.ll2007l ); a 
minority however is composed of detached systems or systems whose eclipses, or at least the 
primary one, show a clear beginning and an end. The small population (~ 60 EBs) around 
V — R ~ 0.5 mag and P ~ 0.3d that is visible in Panel (a), and to a smaller e xtent in Panel 



(b), o f Figure [6] is probably due to foreground contamination. We argue in (jFaccioli et al. 



20071 ) that the period and color of these systems is compatible with them being composed 



of Main Sequence solar mass stars in the Milky Way; the fact that we do not see systems 
with highly eccentric orbit in this foreground population can be attributed to small number 
statistic. 



4.2. LMC results from OGLE-II 

From the 2580 EBs in this sample we selected the 2525 that had a valid V magnitude 
and from these the 2474 EBs for which a phase of secondary minimum was provided; since 
the phase of primary minimum is set to 0, the phase of secondary minimum is equal to our 
4>i _ 4>2] when sec < 0.5 we just took </> sec — > 1 — (p sec . Figure [7| shows the Color Magnitude 
Diagram and Figure [8] shows the Color Period Diagram. The V — Jdia axis interval in both 
diagrams has been chosen so that it has roughly the same range as the V — R axis in Figures. 
|5] and [6] once the difference between V — R and V — I colors is taken into account; because of 
this cut 19 OGLE-II EBs are not shown in Figures [7| and [HI The basic features found in the 
MACHO sample are also evident in the OGLE-II sample, with evolved EBs having mostly 
circular orbits and EBs on the Main Sequence having a broad range of eccentricities. 



4.3. SMC results from MACHO and OGLE-II 



For the SMC the fit, applied to the sample of 1509 EBs introduced in ( jFaccioli et al. 



20071 ). gave acceptable minima determinations in at least one band for 1380 EBs. We observe 
circularization past the Main Sequence in these systems as well as shown by Figures M and 

□3 

From the 1351 EBs in the OGLE-II sample we selected the 1317 that had a valid V 
magnitude and for which a phase of secondary minimum was provided. Figure [IT] shows 
the Color Magnitude Diagram and Figure d2] shows the Color Period Diagram. Again the 
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Fig. 7. — Color-Magnitude Diagram for 2474 EBs in the LMC OGLE-II sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i] < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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Fig. 8. — Color- Period Diagram for 2474 EBs in the LMC OGLE-II sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i| < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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Fig. 9.— Color-Magnitude Diagram for 1380 EBs in SMC MACHO sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i| < 0.6, Panel (c): 0.6 < \<f> 2 - <f>i\ < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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Fig. 10.— Color- Period Diagram for 1380 EBs in the SMC MACHO sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i| < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 



-21 - 



V — Ibia axis interval in both diagrams has been chosen so that it has roughly the same range 
as the V — R axis in Figures. [9] and [10] once the difference between V — R and V — I colors is 
taken into account; because of this cut 15 OGLE-II EBs are not shown in Figures [TT1 and [T2"l 
The basic features found in the MACHO sample are also evident in the OGLE-II sample, 
with evolved EBs having mostly circular orbits and EBs on the Main Sequence having a 
broad range of eccentricities. 



4.4. Period-Phase Difference and Color-Phase Difference Diagrams 

The effects of period on eccentricity are further highlighted in Figure [TBI showing |0i — 02 1 
as a function of period for the four data sets we consider. The LMC EBs from MACHO 
are shown in the upper left panel: the striking dependence of eccentricity on period for 
low periods is clearly shown and two "steps" at about P ~ 0.5d and P ~ 1.5d are clearly 
evident. The lowest period objects in the LMC sample come mostly from contamination by 
foreground sources: there are ~ 60 systems that, in view of their short periods (P < 0.5d) 
and relatively high color ( V — R ~ 0.5 mag) are most likely composed of solar type stars 
located in the Milky Way fjFaccioli et al. 2007 ). This is shown in Table H] that reports data 



for the lowest period objects (defined as P < 0.5d) and shows that the most significant 
contribution to the first "step" at P ~ 0.5d is largely due to this foreground population. 
Figure [IB] suggests that P ~ 1.5d is the likely circularization cutoff period for the LMC 
and P ~ Id is the likely circularization cutoff period for the SMC; these cutoff periods are 
suggested both both by the MACHO samples and the OGLE-II samples. 

The absence of systems with eccentric orbits for P > 20d is due to these systems being 
almost exclusively evolved and therefore having circularized their orbit on the Main Sequence. 
The LMC EBs from OGLE-II are shown in the upper right panel of Figure [13] for the 2474 
EBs for which |0i — 02 1 was provided; these are given to the second decimal place only and 
this accounts for the horizontal "stripes" in the OGLE-II panels in Figures [TBI and [T4l below. 
The features are the same as the MACHO sample: we see again a "step" at P ~ 1.5d where 
there is a sharp increase in the range of eccentricities; EBs with eccentric orbits are almost 
absent for P > 20d; the almost complete absence of EB with P ~ 0.5d probably indicates 
an absence of foreground contamination in the OGLE-II sample. The analogous result for 
the SMC are shown in the lower left panel of Figure [TB] for the MACHO sample and in 
the lower right panel of Figure (TBJ for the OGLE-II sample. The most notable feature of 
these diagrams is the different circularization cutoff: there is a sharp increase in the range of 
eccentricities at P ~ Id as opposed to P ~ 1.5d for the LMC. We explain this difference with 
the fact that EBs in the SMC are younger, on average, than the LMC EBs; therefore only the 
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Fig. 11. — Color-Magnitude Diagram for 1317 EBs in the SMC OGLE-II sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i] < 0.6, Panel (c): 0.6 < \<f> 2 - 0i| < 0.65, Panel 
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Fig. 12.— Color-Period Diagram for 1317 EBs in the SMC OGLE-II sample. Panel (a): 
102 - 0i | < 0.51, Panel (b): 0.51 < |0 2 - 0i] < 0.6, Panel (c): 0.6 < |0 2 - 0i| < 0.65, Panel 
(d): |0 2 -0i| >0.65. 
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shorter period systems, whose components are closer, have had enough time to circularize 
their orbit. The fact that EBs in the two Clouds have, on average, populations at different 
stages of stellar evolution is shown in Table [3J the fraction of EBs belonging to the Main 
Sequence is considerably higher in the SMC than in the LMC. There is also a population of 
~ 9 EBs with 20d < P < lOOd and V — R ~ 0.3 mag — 0.5 mag and strongly eccentric orbits. 
The same qualitative features are shown by the OGLE-II sample; in particular the cutoff 
at P ~ Id is again clearly visible. The OGLE-II SMC sample also reveal a sizeable (~ 20) 
population of long period P > 80d, evolved (1 mag < V — Jdia ^ 2 mag) objects which are 
absent from the MACHO sample; several of these systems have eccentric orbit; their light 
curves reveal that many of them are detached systems. We attribute the fact that we did 
not find as many long period objects in the SMC MACHO sample to the fact that MACHO 
observed the SMC less frequently than the LMC, thus making less likely to find long period 
objects. 

The effects of color on eccentricity are shown in Figure [141 showing |0i — 02 1 as a function 
of V-R for the MACHO samples and of V - J D ia for the OGLE-II samples. The LMC EBs 
from MACHO are shown in the upper left panel; the Main Sequence corresponds there 
to —0.2 mag < V — R < 0.2 mag; the LMC EBs from OGLE-II are shown in the upper 
right panel for the 2474 EBs for which \<f>i — <p2 | was provided; in both cases it is evident 
that the bluest systems have the broadest range of eccentricities. The SMC EBs from 
MACHO are shown in the lower left panel; the Main Sequence there corresponds again; to 
—0.2 mag < V — R < 0.2 mag the LMC EBs from OGLE-II are shown in the lower right 
panel: in both cases is evident that the highest eccentricities are found in the bluest systems. 



5. Discussion 

We assessed the significance of our results, in both Clouds, both for the MACHO and 
the OGLE-II samples. In each case we proceeded in the following way. First we selected 
two subsamples containing EBs with circular orbits and EBs with highly ecce ntric orbits 



respe ctively, and compared their colors via a Kolmogorov-Smirnov (KS) test (IPress et al. 



19921 ). This test estimates the probability for two sets of random values to be drawn from 
the same distribution. In our case we compare the distributions of colors, which we take 
as a rough proxy for evolutionary status since redder systems are generally more evolved. 
This allows us to estimate the probability that, on average, the two subsamples contain 
systems in the same state of stellar evolution. Figure [13] shows that the cutoff period above 
which elliptical orbits start to appear is P ~ 1.5d for the LMC and P ~ Id for the SMC. 
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Fig. 13. — Period-Phase Difference Diagram for the four datasets considered. 
Upper left: LMC MACHO sample (4510 EBs). 
Upper right: LMC OGLE-II sample (2474 EBs). 
Lower left: SMC MACHO sample (1380 EBs). 
Lower right: SMC OGLE-II sample (1317 EBs). 

For the LMC the circularization cutoff is evident at P ~ 1.5d; for the SMC it is evident at 
P ~ Id. 
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Fig. 14. — Color-Phase Difference Diagram for the four datasets considered. 
Upper left: LMC MACHO sample (4510 EBs). 
Upper right: LMC OGLE-II sample (2474 EBs). 
Lower left: SMC MACHO sample (1380 EBs). 
Lower right: SMC OGLE-II sample (1317 EBs). 
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Furthemore above P ~ 20d there are again almost only EBs with circular orbits in both 
Clouds. Therefore we concentrated on the range 1.5d < P < 20d for the LMC and Id < P < 
20d for the SMC. In these ranges we selected both EBs with circular orbits (|0i — ^2 1 < 0.51), 
and EBs with eccentric orbits (|0i — 02 1 > 0.6 =>- e > 0.16). The cutoff of 0.51 for circular 
orbits was chosen because, as Figures [3] and H] show, the error in the phase difference is on 
average 0.01 or less; thus a cutoff at 0.51 usefully discriminates between circular and eccentric 
orbits; for the OGLE-II samples the cutoff is also appropriate since the phase of secondary 
minimum is reported up to two decimal places. For the MACHO samples we used the V — R 
color; for the OGLE-II samples we used both the V — I color and to the V — Jdia color. 

For the purpose of this discussion we call these subsamples the all ages subsamples 
because they contain EBs of all ages; in particular they contain EBs that have evolved past 
the Main Sequence. We then selected two more subsamples of blue and young EBs, again 
comprised of EBs with circular orbits and EBs with highly eccentric orbits. These subsamples 
were selected with a color cut whose exact definition depends on the EB sample (see below), 
and with the same cuts in |0i — <fo| an d P as the all ages subsamples defined above and 
their colors were again compared via a KS test. We call these subsamples the young EB 
subsamples because the color cut ensures that they contain only young and unevolved stars 
that only recently settled on the Main Sequence. Table [5] summarizes the properties of these 
subsamples; as t he numbers show all subsamples are large enough for the KS test to be 



safely employed (IPress et eil.l Il992f ) : the values of the KS statistic D for all subsamples are 
reported in Table [6J 

The data reported in Table [6] allow us to draw the main conclusion of this work, which 
hold for both Clouds, namely: 

1. There is a significant probability for the two young EB subsamples of being drawn 
from the same color distribution. 

2. The probability for the all ages subsamples of being drawn from the same distribution 
of colors is vanishingly small. 

These two findings allow us to state that, in both Clouds, EBs start their lives with a broad 
distribution of orbital eccentricities and circula rize their orbits as their evolve on and past the 



Main Sequence. This result is consistent with ( jLucy fc Sweeney! 1 19711 ) who find that, in their 



sample of 103 spectroscopic binaries, long per i od sy stems with late type giant components 



have mostly circular orbits; iLucy fc Sweeney! (1197ll ) suggest that tidal interactions due to 



expansion after the Main Sequence phase of one or both components is responsible for their 
circular orbit. We note finally that the use of the V — I cut in the OGLE-II samples gives 
results slightly more consistent with the results we obtain from the MACHO samples than 
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the use of the V — /dia cut, although our conclusion does not qualitatively change. For 
this reason in the histograms of Figures [15] and [16] we show the V — I color rather than the 
V — Jdia color. 



5.1. Influence of a third light 

One concern with our procedure is the possible influence of a third light. If one EB 
system is close to a bright and blue undetected third star, the overall color of the system is 
skewed toward the blue. This could result in an EB being erroneously selected as "young" 
by the color cuts described above, whereas the system is actually in a more advanced stage 
of evolution and biasing our conclusions. 

We assessed the possible significance of this effect by fitting the V light curves of the 
MACHO LMC EBs selected with the V-R < -0.06 mag cut, with with and without a third 
light as a free fit parameter; we then did the same for the SMC with the V — R < —0.05 mag 
cut. If the V band light curve fit improves considerably by using a third light as a free fit 
parameter, this means that an undetected blue third star is present; such system should 
therefore be excluded from the KS tests because the binary system is actually redder and 
probably more evolved than its V — R suggests. We point out that the only aim of these fits 
is to assess the significance of a possible undetected blue third star, and not to determine 
any physical parameters for the EBs. 

Since the systems we are dealing with are unevol ved we exp e ct them to be detache d and 

We 



therefore to be well described by the EBOP model ([ 



actually performed the fits using the JKTEBOF0 (ISouthworth. Maxted fc Smalleyl 



tzellll98ll: IPopper & Etzel 



198 



2004 



Southworth et al.l 120041 ) code, based on EBOP and which adds several modifications and 
extensions that make it easier to use, especially when fitting a large number of light curves 
We o btained preliminary starting values for the fit parameters using the DEBilJ^I code (D 



;evor 



20051 ). a fitting code suited to automatically analyze large numbers of detached EBs. We 
then refined the fits using JKTEBOP; this was done by performing a first fit using the 
parameters computed by DEBiL as starting values, then eliminating those points at least 3 
standard deviations away from this first fit and redoing the fit; this second step allowed in 
many cases to considerably improve the quality of the fit by excluding outlying points. We 
used the following criterion to quantify whether or not the quality of the fit improved as a 



http : //www . astro . keele . ac . uk/~{} j kt/ codes/ j ktebop . html 
http : //www . cf a . harvard . edu/~{} j devor/DEBiL . html 
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Table 4. Summary of very low period objects in the MACHO LMC sample. 



Population 


Number of objects 


Very low period objects a 


83 


Very low period objects' 1 with circular orbit b 


56 


Foreground population 


62 


Very low period objects a in the foreground population 


40 


Very low period objects 8, with circular orbit b in the foreground population 


36 



a Defined as P < 0.5d. 
b Defined as |0i - 2 | < 0.51. 



Table 5. Summary of MACHO and OGLE-II EB eccentricity data. 



Galaxy 


Sample 


Circular' 


1 Eccentric b 


Circular Young c 


Eccentric Young c 


Period Range 


LMC 


MACHO 


1537 


237 


260 


79 


1.5d < P < 20d 


LMC 


OGLE-II (/ mags.) 


1084 


119 


145 


30 


1.5d < P < 20d 


LMC 


OGLE-II (Jdia mags.) 


1084 


119 


211 


44 


1.5d < P < 20d 


SMC 


MACHO 


724 


79 


104 


20 


Id < P < 20d 


SMC 


OGLE-II (I mags.) 


468 


70 


212 


39 


Id < P < 20d 


SMC 


OGLE-II (Jdia mags.) 


468 


70 


233 


45 


Id < P < 20d 



a Defined as |</>i - <f> 2 \ < 0.51. 
b Defmed as \(j>i -<j> 2 \> 0.6. 



c Defined as V-R < -0.06 mag for the LMC MACHO sample, as V-R < -0.05 mag for the SMC MACHO 
sample, and as V — I < mag or V — Jdia < mag for both OGLE-II samples. 
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Fig. 15. — Color distribution for EBs with no detected eccentricity (continuous line) and 
for EBs with eccentricity > 0.16 (dashed line) with 1.5d < P < 20d for the LMC and 
Id < P < 20d for the SMC. Upper left: LMC EBs from MACHO. Upper right: LMC EBs 
from OGLE-II. Lower left: SMC EBs from MACHO. Lower right: SMC EBs from OGLE-II. 
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Fig. 16. — Color distribution for young EBs with no detected eccentricity (continuous line) 
and for EBs with eccentricity > 0.16 (dashed line) with 1.5d < P < 20d for the LMC 
and Id < P < 20d for the SMC. Upper left: LMC EBs from MACHO: V-R < -0.06 mag. 
Upper right: LMC EBs from OGLE-II: V-I < mag. Lower left: SMC EBs from MACHO: 
V-R < -0.05 mag. Lower right: SMC EBs from OGLE-II: V - I < mag. 
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result of fitting for a third light. The typical number of points in the blue MACHO light 
curves is ~ 690 for the LMC and ~ 840 for the SMC and the number of fit parameters 
is ~ 10; we therefore assumed a typical number of degrees of freedom (dof) for our fits of 
800. For 800 dof, the probability of getting X 2 > 960 x 2 /dof > 1.2 by chance is < 10~ 4 
whereas the probability of getting x 2 > 800 =>- x 2 /dof > 1 by chance is 0.4933. Therefore, 
if for a V light curve we have x 2 /dof > 1.2 without fitting for a third light and x 2 /dof < 1 
by fitting for it, we conclude that a blue third object is present; the EBs is then excluded by 
the young EB subsamples and the KS tests are performed again on them. This procedure 
gave the subsamples described in Table the results of the KS tests are reported in Table 

E 

As Table [5] shows, the probability of the two subsamples of being drawn from the same 
color distribution remain significant after accounting for the possible influence of a third 
light. 



5.2. Results from EROS 



We co nsidered the sample of 79 EBs in the bar of the LMC found by the EROS col- 
laboration (IGrison et al.lll995l ). We cross correlated this sample with both the MACHO 
sample, finding 42 matches, and with the OGLE-II sample finding 54 matches; therefore this 
sample overlaps only partly with the two larger samples and we performed the same analysis 
described above, finding comparable results, despite much lower number statistic. Of the 79 
EBs in the sample we selected, in the 1.5d < P < 20d period range, 33 EBs with circular 
orbits and 9 EBs with eccentric orbits. We then applied the KS test to the distributions of 
the B magnitudes for t hese two samples, since, despite their small size, the test can still be 
used ( IPress et al.lll992l ). We found that the two distributions are different at ~ 78% level; 
the fact that the confidence level is not higher can be attributed to small number statistic. 
We did not attempt to select subsamples of young EBs since their numbers would have been 
too small for the KS test to be applied. 



6. Conclusion 

We have presented a study of orbital circularization in the Magellanic Clouds using EBs 
samples compiled both by the MACHO and by the OGLE-II collaboration, as well as a sample 
for the bar of the LMC compiled by the EROS collaboration. We have shown that in the 
LMC binary stars with period in the range 1.5d— 20d with initially eccentric orbits circularize 
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Table 6. Results of KS tests for the subsamples of Table [U 



Galaxy 


Sample 


Subsamplc 


KS Statistic D 


Probability of D 


LMC 


MACHO 


All ages 


0.256 


2.3 x 10~ 12 


LMC 


OGLE-II 


All ages (I mags.) 


0.351 


3.3 x 10~ 12 


LMC 


OGLE-II 


All ages (7dia mags.) 


0.362 


5.7 x 10~ 13 


LMC 


MACHO 


Young 


0.092 


0.661 


LMC 


OGLE-II 


Young (/ mags.) 


0.131 


0.757 


LMC 


OGLE-II 


Young (7 DIA mags.) 


0.157 


0.304 


SMC 


MACHO 


All ages 


0.216 


0.002 


SMC 


OGLE-II 


All ages (J mags.) 


0.141 


0.162 


SMC 


OGLE-II 


All ages (7 D ia mags.) 


0.219 


0.005 


SMC 


MACHO 


Young 


0.156 


0.776 


SMC 


OGLE-II 


Young (/ mags.) 


0.159 


0.354 


SMC 


OGLE-II 


Young (J D ia mags.) 


0.186 


0.121 



Table 7. Summary of MACHO eccentricity data for EBs with negligible influence of a 

blue third light. 



Galaxy 


Sample 


Circular a Young c 


Eccentric b Young c 


Period Range 


LMC 


MACHO 


240 


78 


1.5d < P < 20d 


SMC 


MACHO 


99 


19 


Id < P < 20d 



a Dcfmed as \cj>i - <j> 2 \ < 0.51. 
b Dcfmed as |0i -<£ 2 | > 0.6. 

c Defined as V-R < -0.06 mag for the LMC MACHO sample, and as V-R < 
-0.05 mag for the SMC MACHO sample. 



Table 8. Results of KS tests for the subsamples of Table [71 



Galaxy Sample Subsample KS Statistic D Probability of D 



LMC MACHO Young 0.098 0.597 

SMC MACHO Young 0.153 0.816 
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during post Main Sequence evolution, with a cutoff in color at about V — R < 0.2 mag: our 
data are consistent with stars starting their lives with a wide range of eccentricities and 
circularizing their orbit as they evolve, both on the Main Sequence and in their red giant 
phase. For the SMC the same conclusion holds, for a period range Id— 20d. We have obtained 
the same results with independently assembled datasets, thus enhancing our confidence in 
the validity of our conclusions. 



This work uses public domain data from the MACHO project whose work was performed 
under the joint auspices of the U.S. Department of Energy, National Nuclear Security Ad- 
ministration by the University of California, Lawrence Livermore National Laboratory under 
contract No. W-7405-Eng-48, the National Science Foundation through the Center for Parti- 
cle Astrophysics of the University of California under cooperative agreement AST-8809616, 
and the Mount Stromlo and Siding Spring Observatory, part of the Australian National 
University. KHC's work is performed under the auspices of the U.S. Department of Energy 
by Lawrence Livermore National Laboratory in part under Contract W-7405-Eng-48 and in 
part under Contract DE-AC52-07NA27344. This work uses public domain data obtained 
by the OGLE project. LF acknowledges the kind hospitality of the Institute of Geophysics 
and Planetary Physics at Lawrence Livermore National Laboratory and of the Harvard- 
Smithsonian Center for Astrophysics where part of the work was done. We also thank the 
referee for helpful sugge stions, and we thank Leon Lucy for pointing out to us the reference 
(ILucy &: Sweeney! 1 197ll ). 
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